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ABSTRACT 

We reassess constraints on the cosmological baryon density from observations of the mean 
decrement and power spectrum of the Lyman-a forest, taking into account uncertainties in all 
free parameters in the simplest gravitational instability model. The uncertainty is dominated 
by that of the photoionizing background, but incomplete knowledge of the thermal state of the 
intergalactic medium also contributes significantly to the error-budget. While current estimates 
of the baryon fraction from the forest do prefer values that are somewhat higher than the big bang 
nucleosynthesis value of Vl^ti 2 = 0.02 ± 0.001, the discrepancy is at best about 3a. For instance, 
assuming the highest estimate of the ionizing background, as indicated by recent measurements 
of a large escape fraction from Lyman-break galaxies by Steidel, Pettini & Adelberger, we find 
ilbh 2 — 0.045±0.008. A recent measurement of the ionizing background from the proximity effect 
by Scott et al., on the other hand, implies Q,bh 2 = 0.03 ± 0.01. We provide an expression from 
which future likelihoods for flbh 2 can be derived as measurements of the ionizing background 
improve - consistency among constraints from the forest, nucleosynthesis and the microwave 
background will provide a powerful test of the gravitational instability model for the forest, and 
for large scale structure in general. We also develop a formalism which treats lower bounds on the 
baryon density in a statistical manner, which is appropriate if only a lower bound on the ionizing 
background is known. Finally, we discuss the implications of the escape fraction measurement 
for the age, structure and stellar content of Lyman-break galaxies. 



Subject headings: cosmology: theory - intergalactic medium - large scale structure of universe; 
quasars - absorption lines 



1. Introduction 

It has long been recognized that the intergalactic medium (IGM) is highly ionized, and likely contains a 
substantial fraction of the universe's baryons at redshifts z ~ 3 (e.g. Gunn & Peterson 1965, Steidel & Sargent 
1987, Fukugita, Hogan & Peebles 1998 and references therein). Constraints on the baryon density f^ft 2 
(where fit is a fraction of the critical density, and h is the Hubble constant today in units of 100 km/s/Mpc) 
from observations of the Lyman-a (Lya) forest rely on assumptions about the size and geometry of the 
absorbing structures (see e.g. Rauch & Haehnelt 1995). Recent advance in theoretical modeling of the forest 
provides a framework in which fluctuations seen in absorption arise naturally out of gravitational instability, 



and which successfully matches the gross observed properties of the forest (e.g. Bi, Borner & Chu 1992; Ccn 
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ct al. 1994; Zhang et al. 1995: 


Reiscncggcr & Miralda-Escudc 199E; Hcrnquist et al. 1996; Miralda-Escudc 


ct al. 1996|; Mucckct et al. 1996 


; 


Bi & Davidsen 1997; 


Bond & Wadsley 1997: 


Croft ct al. 1998; Hui, Gncdin 


fe Zhang 1997; Thuens et al. 1999; Bryan et al. 1999; Croft et al. 1999; Nusser & Hachnelt 200C; but see also 



Meiksin, Bryan & Machacek 2001 for possible problems, and Zaldarriaga et al. 2001 for consistency checks ). 
This framework makes definite predictions about the spatial distribution of the absorbing structures, making 



it possible to obtain useful constraints on fli,h 2 from forest observations (e.g. Rauch ct al. 1997; Weinberg et 
al. 1997| ; [McDonald et al. 200C| ) - once the photoionizing background, which determines the overall neutral 
fraction, is specified. 

Three recent developments motivate us to re-examine the Lya forest constraints on the cosmological 
baryon fraction. First, the recent microwave background anisotropy measurements favor a baryon fraction 
that is higher than the big bang nucleosynthesis value inferred from deuterium measurements (the former 
gives 0.022 < flt,h 2 < 0.040, 2-er limits for a fiat scale-invariant inflation model, see de Bernardis et al. 2000]; 



Hanany, S. et al. 2000; Tegmark fc Zaldarriaga 2000; Jaffe et al. 2000; Lange et al. 2001; the latter gives 



0.02±0.001, 1-cr error, see Buries, Nollet & Turner 2000; see also Buries & Tytler 1998 for earlier constraints). 
Second, Steidel, Pettini & Adelberger (2001) reported measurements of a large escape fraction of ionizing 
photons from Lyman break galaxies (LBG's) at z ~ 3, which imply an extragalactic ionizing flux that is 
significantly larger than what one might expect from quasars alone. A larger photo-ionizing background 
implies a smaller neutral fraction, and so a larger baryon density is required to match the given observed 
mean absorption seen in the forest. Lastly, analysis of the Lya forest has moved beyond comparisons with a 
restricted set of cosmological models to a point where several parameters specifying the mass power spectrum 
and the thermal state of the IGM can be determined at the same time from the data (e.g. Choudhury, 
Srianand & Padmanabhan 2000, Zaldarriaga, Hui & Tegmark 2000). Since earlier forest constraints on 
the baryon fraction focused on fixed cosmological models and/or reionization history (and therefore fixed 
thermal state) it is timely to examine how uncertainties in these parameters impact the estimates for 
fit, ft 2 . Because the Lya constraint on the baryon fraction is completely independent of constraints from 
nucleosynthesis or the microwave background, comparing all three provides a powerful consistency check of 
the inflation + cold dark matter structure formation model, as well as a test of current ideas about the 
nature of the forest. 



We briefly review in §2.1 the gravitational instability model for the Lya forest, and how it is used to 
obtain estimates of the baryon density. Discussions in the literature have focused on two different approaches. 
One (e.g. Haehnelt et al. 2000) implicitly assumes that most of the universe's baryons is in the IGM 0, and 
that the total ionizing background can be largely accounted for by directly summing up the contributions 
from known sources, or obtained from indirect measurements of the sum total background, such as those 
using the proximity effect. The other approach (e.g. Rauch et al. 1997; Weinberg et al. 1997) is more 
conservative, assuming only a lower limit to the ionizing background from known sources - this gives a lower 
bound on the neutral fraction and therefore a lower bound for the f2f,ft 2 . We will discuss both approaches. 



We describe in §2.2.1 how lower bounds given in the literature can be interpreted in a statistical manner, 
generalizing earlier work where only strict lower bounds are discussed. As we will see, uncertainties in the 



A notable exception is Weinberg et al. (1997) who attempted to give the most conservative lower bounds on the baryon 
density. We provide here a statistical framework for interpreting such lower bounds. 

2 Estimates so far do give values for the baryon density that are close to or even higher the the nucleosynthesis value (e.g. 
Rauch et al. 1997, McDonald et al 2000). Detailed accounting of the universe's baryon budget also supports the notion that 
most of the baryons are in the IGM (Fukugita et al. 1998). 
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level of the ionizing background dominate uncertainties in the current bound on f^/i 2 . In §2.2.2, we relate 



the lower-bound likelihood in §2.2.1 to the (differential) likelihood for baryon density, assuming effectively 
that most of the universe's baryons are in the IGM and that the ionization background is known. As we will 
emphasize in §2.2.2, non- negligible uncertainties in the estimates for f2f,/i 2 remain even if the ionizing back- 



ground is known to high accuracy. This is because of degeneracy with other parameters in the gravitational 
instability model. We will examine which parameters are the most important in this regard. The likelihood 
of Clbh 2 for a fixed ionizing background is also given (an explicit expression is also presented in §|J), which 
can be used to compute future likelihoods as measurements of the ionizing background improve. 

Some of the conclusions in the above sections depend critically on the ionizing flux inferred from ob- 
servations of Lyman break galaxies (Steidel et al. 2001). The fraction of hydrogen-ionizing photons that 
escape from these galaxies is significantly higher than what is observed in local galaxies, and what is naively 
expected in simple theoretical models. In §^, we investigate the type of cold gas distribution that would 
allow these early galaxies to have a high escape probability for ionizing photons. Finally, we conclude in §||, 
where the issue of escape fraction from quasars is also discussed. 

Some of the issues in this paper have been addressed in two recent works. Haehnelt et al. (2000) used 
previous bounds to derive new limits on flbh 2 , using the Steidel et al. measurement of the escape fraction 
- the treatment assumed the same cosmological model and reionization history upon which the previous 
respective bounds were based. They also briefly discussed the significance the Steidel et al. results in the 
context of population synthesis model spectra, and the escape fraction expected based on the amount of gas 
inferred from the star-formation rate. Choudhury, Srianand & Padmanabhan (2001) used the lognormal 
approximation to derive bounds on the ratio of (flbh 2 ) 2 to the ionizing flux, where the thermal state of the 
gas is allowed to vary, for a fixed cosmological model. Here, we examine the f^/i 2 constraints using N-body 
simulations for a suite of cosmological models (i.e. varying power spectrum) and thermal histories. We also 
consider the implication of the high observed escape fraction for the age, structure and stellar content of 
Lyman-break galaxies in greater detail. 



2. Baryon Fraction from the Lyman a Forest 

2.1. The Gravitational Instability Model for the Forest 

We describe briefly here the Lya forest model used in this paper. It is largely motivated by the success of 
numerical simulations in matching the observed properties of quasar absorption spectra (see Cen et al. 1994 
and related references in §0). The Lya optical depth r along the line of sight is given by A(l + 5) a , mapped 
to redshift space, and smoothed with the thermal broadening window (see e.g. Hui, Gnedin & Zhang 1997 
for details). The quantity S is the baryon overdensity Sp/p. The constant a is given by 2 — 0.7(7 ~ 1) where 
7 is the equation of state index - i.e. the temperature of the gas T is related to density by T = Tq(1 + <5) 7_1 , 
where To is the temperature at mean density (Croft et al. 1997, Hui & Gnedin 1997). The constant A is 
given by (e.g. Weinberg et al. 1997) 3 
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3 A = nmcH(z) 1 cr ct o, where nni is the neutral hydrogen density at 5 = 0, c is the speed of light, and H(z) is the Hubble 
parameter at redshift z, and <r a o = 4.478 X 10~ 18 cm 2 : the Lya absorption cross-section is given by cr a Q times the Voigt profile. 
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0.76 



(188 



4.461 



y/Q m (l + z) 3 + n A + (1 - Q m - 00(1 + z) 5 



where /i is the Hubble constant today in units of 100 km/s/Mpc, z is the redshift, Xh is the mass fraction 
of baryons that is in hydrogen, X^ c is the same for helium, the assumption being that both H and He are 
highly ionized. [] The factor of 4.461 in the last term corresponds to a model with Q m = 0.3, Q\ = 0.7 and 
z = 3. We have used the recombination rate coefficient a roc . = 4.2 x 10~ 13 (T/10 4 K)~°- 7 cm 3 s _1 (Rauch et 
al. 1997). 



The dimensionless quantity Jhi is defined by 



Jhi x 10 21 ergs 1 cm 2 Hz 1 ster 1 
/ 47rj 1 ,cr H i(^) — / / 47rcr H l(^) — 



(2) 



/3 + 3 



where j v is the specific intensity of the ionizing background, i/jji = c/912A, ohi is the ionizing cross-section 
for HI (<7hi(V) is approximately 6.3 x 10 _18 cm 2 (v/vm)~ 3 for v just above i/hi; but see below), and f3 is the 
slope of j v just above ^hi: ii/ oc j/ - ' 3 . The photoionization rate Thi is related to Jhi by 



r H i = 4.3 x 10~ 12 J H iS 



(3) 



The expression here, together with eq. (g) above, gives the correct relationship between j„ HI and Thi that 
is accurate to within 3%, for all values of (3 between and 3. This takes into account slight departure of 
chi(^) from an exact v~ 3 power-law (see Osterbrock 1989, eq. 2.4). 

Given a structure formation model such as the Cold Dark Matter (CDM) model, the distribution of 
6 follows from gravitational instability, and can be predicted using hydrodynamic, N-body (with suitable 
smoothing) or Hydro-PM (Gnedin & Hui 1998) simulations. Combining the 8 fluctuations with the relation 
for the optical depth r = A(l + 5) a (with suitable mapping to redshift space and thermal broadening), 
one can predict a host of observable properties for the Lya forest, for any given value of A and a. Two 
observables that have attracted a lot of attention and have been well measured are the transmission power 
spectrum and the mean decrement (1 — (e~ T )). These measurements can be used to yield constrains on 
parameters in the model, including A which is proportional to (fl^h 2 ) 2 . 

In this paper, we make use of the analysis by Zaldarriaga, Hui & Tegmark (2000), who obtained 
constraints on 6 different parameters by matching the observed transmission power spectrum and mean 
decrement at z — 3 from McDonald et al. (2000). The 6 parameters are the mass power spectrum slope and 
normalization, the gas temperature at mean density To, the equation of state index 7, the normalization 
factor A (eq. [Q), and fc/, which is the smoothing scale that defines the smoothing of gas relative to dark 
matter. f\ The quantity of interest for our purpose here is the combination AT® 7 = A(Tq/10 4 K ) 7 , which 
is directly proportional to (Vl^ti 2 ) 2 / Jhi- From the analysis of Zaldarriaga et al. (2000), it is straightforward 
to obtain the likelihood for AT®- 7 ', marginalized over the other parameters. This is shown in Fig. 0a. It 



4 Assuming instead that He is only singly ionized, which amounts to changing 2Xn c to Xn c , will make little difference given 
our other sources of uncertainties. Whether helium is singly ionized or not of course has profound implications for the spectrum 
of the ionizing radiation — that is taken into account by the parameter Jhi- 

5 In principle, parameters such as Q m and also affects the predicted forest properties, but at redshift 3, Oa is negligible 
and Q m is close to unity. The cosmological constant does affect the translation between velocities and distances - the mass 
power spectrum slope and normalization are essentially fixed at the relevant velocity scales. 
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provides the starting point for our analysis below. It is worth mentioning that the mean decrement from 
McDonald et al. (2000) has not been corrected for a possible bias introduced by continuum-fitting. However, 
comparing their measurement with that from Rauch et al. (1997), who did attempt a correction (which is 
model dependent), reveals that the bias is comparable to the error-bar given by McDonald et al., who gives 
(e- T ) = 0.684 ± 0.023 at z = 3. 



2.2. Likelihood for the Baryon Fraction 

The likelihood for AT®' 1 shown in Fig. [I]a can, in principle, be translated directly into a likelihood for 
£lbh 2 after marginalization over other parameters that appear in A (eq. jjj), including Jhi, h, etc. However, 
as mentioned in §[l| there are at least two reasons why, strictly speaking, one obtains only lower bounds for 
the baryon fraction - first, not all baryons are necessarily in the forest, although a large fraction does appear 
to be (i.e. fl^h 2 in eq. (IJ refers to only those baryons in the forest); second, one often has strictly speaking 
only lower limits on the ionizing background Jhi, based on summing over contributions from known sources 
(an exception is the use of the proximity effect; another exception is an observational upper limit on Jhi 
from searches for fluorescent Lya emission; both are discussed below). In this section, we will discuss both 
the likelihood for flbh 2 - as if all or most baryons arc in the forest, and the total Jhi is known - as well as 
the lower-bound- likelihood for flbh 2 . 



2.2.1. A Probabilistic Analysis of Lower Bounds 

Let us begin by developing the idea of a lower-bound-likelihood. Previous work often gave lower bounds 
on the baryon fraction as if they were strict bounds. Of course, uncertainties in Jhi, in the measured mean 
decrement, and in other parameters such as To imply that such bounds have to be given a probabilistic 
interpretation. We show here the proper quantity to consider is a lower limit on the probability that flbh 2 is 
larger than some value. 

We have from eq. (0) Clbh 2 = ^/ ' AT^' 1 Jml '(<?) , where q represents the parameters h, f2 m , Cl\, and f(q) 
is some function which can be read off from eq. ([j]) (we will focus on z = 3 in this paper). The probability 
that flbh 2 is larger than some value Bq is given by 

i 

P(n b h 2 )d(Q b h 2 ) (4) 
= / dqd(AT^)P(q)P(AT°- 7 ) f P,u(Jm)dJ m 
> f dqd(AT°- 7 )P(q)P(AT°- 7 ) f P(Jm)dJ m 

= Pyiilbh 2 > B ) 



B„ 



where P(q), P(AT® - 7 ), P{ Jhi) represent the probability distributions of q, AT®- 7 and Jhi- We use P a ii(Jffl) 
to denote the probability distribution for the ionizing intensity from all possible sources, while the probability 
distribution P( Jhi) that we will use, at least for some of the calculations below, comes from a summation 
over contributions from known sources, including Lyman-break galaxies and quasars. The inequality above 
can be understood as follows. For any random variable Z = X + Y with X, Y > and independent, we have 
J~Pz(Z)dZ = J ° c dYP Y (Y)£_ Y dXP x (X) > J cc dYP Y (Y)J^dXP x (X) = J™dXP x (X). By the 
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same reasoning, the fact that the amount of baryons in the forest is necessarily smaller than the cosmological 
total only serves to strengthen the inequality in eq. We introduce the symbol P>(f2b/i 2 > Bq) to denote 
the lower limit to the cumulative probability that fl^h 2 is larger than some value B . 

There are therefore, three probability distributions that we have to specify: P(AT® ) is already given in 
Fig. ||a - that is the input from observations of the Lya forest; P(q) does not actually affect our conclusions 
very much, but we will assume a flat universe, with independent Gaussians for £l rn h (0.28 ± 0.09; Eisenstein 
& Zaldarriaga 2001) and h (0.72 ± 0.08; Freedman et al. 2001) (at z = 3, Oa plays essentially no role); 
-P(Jhi) is what we turn to next. There are three main sources of information. 

A recent analysis of proximity effect by Scott et al. (2000) (see also Giallongo et al. 1996, Cooke, 
Epsey & Carswell 1997) gave Thi = 1.9±i o x 10- 12 s-\ averaged from z = 1.7 to 3.8, which translates into 
Jhi = 0.44^Q 23 s _1 . This is presumably the ionizing background from all possible sources. 

Steidel et al. (2001) gave j Vm — 1.2 ± 0.3 x 10~ 21 erg s _1 cm -2 Hz -1 sr _1 from the Lyman-break galaxies 
(LBG's) at z ~ 3. | This is obtained by examining the spectra of those LBG's that are in the bluest quartilc 
of observed LBG UV colors. As pointed out by Steidel et al. (2001), a conservative interpretation would 
therefore be that no ionizing radiation escapes from the other galaxies, in which case j Um would be a factor of 
4 smaller: j Um = 0.3±0.075 x 10 _21 ergs~ 1 cm~ 2 Hz -1 sr _1 . There are in addition a number of uncertainties, 
including uncertainty in the mean absorption by the intergalactic medium, and the use of the luminosity 
function down to magnitudes that are fainter than those used for the escape fraction measurements (see 
Steidel et al. for details). These uncertainties are smaller than the factor of 4 uncertainty in interpretation. 
To convert the above into Jhi, we use = 1.4 (see eq. [|; we will discuss the reason for this choice of in 
fj- 

Finally, there is the contribution from quasars. Rauch et al. (1997) discussed the various measurements 
and calculations in some detail, and they arrived at a conservative lower limit of Thi > 7 x 10~ 13 s _1 
(Jhi > 0.16) at z = 2 — 3, when they consider only known quasars with no extrapolation to the faint end of 
the luminosity function. We will treat this as a strict lower limit for Jhi- 

We will work out the lower-bound likelihood in eq. (||) for three different P( Jhi)'s: 



= 1.2 ± 0.3 x 10~ 21 ) from LBG's and Rauch et al.'s 
0.3 ± 0.075 x 10~ 21 ) from LBG's and Rauch et al.'s 

The corresponding quoted error-bars are used in setting the probability distributions for Jhi. In all cases, 
quantities that are by definition positive, e.g. Jhi, Q m h, etc, are constrained to be so. 



6 The quoted number was computed by Steidel et al. assuming Einstein de-Sitter cosmology, but the value is in fact strictly 
independent of cosmology. 

7 In other words, the probability distribution of Jjji from all quasars P a n q (Jhi) is assumed to satisfy the inequality 
P 3,n q .{Jm)d.Jm > ^known q. (Jm)dJm, where J° : can take any given value, and P known q . (Jhi) is the probabil- 

J HI J HI 

ity distribution of Jjji from known quasars, and is modeled as a delta function which peaks at Jjji = 0.175. 



• 1. Jhi from the proximity effect; 

• 2. Jhi from summing the high value of j„ HI (J, 
value from quasars; 

• 3. Jhi from summing the low value of j„ HI (J UB 
value from quasars. 
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We show in Fig. [j]b the result of combining the likelihoods for AT®' 7 , q and Jhi as described above 
using eq. (Q). To reiterate, the main observational inputs are measurements of the mean decrement and 
transmission power spectrum of the Lya forest at z = 3. The y-axis is the lower limit to the likelihood 
that Q b h 2 is larger than some value B , P > (p, b h 2 > B ) as defined in eq. (^). The current big bang 
nucleosynthesis value inferred from deuterium measurements is VL b h 2 = 0.02 ±0.001 (Buries, Nollett, Turner 
2000). This is ruled out at 79%, 99.8% and 96% (or better) for the choices 1, 2 and 3 above for P(J H i). 

The large difference between the curves in Fig. 0b illustrates that the dominant uncertainty by far is 
that due to Jhi- We emphasize, however, even assuming the highest Jhi (model 2), the nucleosynthesis value 
is not strongly excluded by Lya forest measurements, although there is an inconsistency at the ~ 3er level. 



2.2.2. Further Investigations of the Baryon Likelihood 

It is instructive to show a different version of Fig. |l|b, namely the derivatives of those same curves 
(with a negative sign) i.e. —dP > (£l b h 2 > Bq)/oIBq. This is done in Fig. ||a, whose y-axis now represents 
the differential probability distributions for il b h 2 {P(VL b h 2 ), rather than P>{£l b h 2 > Bq); see eq. Q), if one 
assumes that the majority of the baryons is in the forest, and that the probability distributions for Jhi used 
above actually approximate the true distributions for the total Jhi- Phrased in this way, the mean values 
with dispersions for Q b h 2 are 0.03 ±0.01, 0.045 ±0.008 and 0.028 ±0.005 respectively for the choices 1,2 and 
3 for P(Jhi) laid out previously. 

The above represents another way of saying that the present measurements have a weak preference for 
Q b h 2 somewhat larger than the nucleosynthesis value. The error-bars at the moment are dominated by un- 
certainties in Jhi- However, it would be useful to derive a likelihood for fl b h 2 assuming Jhi is exactly known. 
Future improvements in measurements of Jhi can then be directly translated into improved distributions for 
fl b h 2 . To be concrete, we would like to compute P»(P ): 

P*(B' ) ee P(n b h 2 = B' , Jhi = 1) = J dq^P(q)P(AT°- 7 = B' 2 /f(q)) (5) 

which is obtained from eq. (^) by setting P(Jhi) = <5(Jm — !)• We show P*(P = fl b h 2 ) in Fig. ||b. It is 
well-approximated by a Gaussian with a mean of 0.046 and a dispersion of 0.0061. This is the likelihood for 
Q b h 2 if Jhi were exactly unity. If Jhi were 0.5 for instance, then the likelihood for fl b h 2 would be a Gaussian 
with a mean of V0.5 x 0.046 = 0.033 and a dispersion of \/0.5 x 0.0061 = 0.0043. More generally, given any 
P(Jhi), the differential probability distribution for Q b h 2 can be easily computed using P*: 

P{Q b h 2 ) = f -!=P*(B' = n b h 2 /y/l^)P{J m )dJ m (6) 
J v Jhi 

Or, more rigorously, in cases where P(Jm) accounts for only contributions from known but not all possible 
sources, and keeping in mind that not all baryons are necessarily in the forest, one can obtain (by integrat- 
ing eq. ||) a lower bound on the cumulative probability that Q b h 2 is larger than some value - i.e. use 



J B P(£l b h 2 )d(Q b h 2 ) to obtain P > (Q b h? > B a ), just as we have done in §2.2.1. 

Fig. |^b illustrates another important point: substantial uncertainties in fl b h 2 remain even if Jhi were 
known at high accuracy. The uncertainties can be traced back to the spread in AT®' 7 , shown in Fig. |l|a. 
Checking each of the parameters that were marginalized in obtaining P(AT®' 7 ), we find that the main 
sources of errors are, in order of importance, To the temperature at mean density, 7 the equation of state 
index, and kf the smoothing scale of the gas distribution - fixing each of these parameters individually 
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would decrease the variance in AT®- 7 by 43%, 39% and 31% respectively. Interestingly, fixing the power 
spectrum normalization and/or slope, or the mean transmission value (1— mean decrement) decreases the 
variance of AT® 7 by no more than 11%. The latter in particular implies that the current accuracy of the 
mean transmission/decrement measurements at z ~ 3 is not the main limiting factor for the the accuracy 
of the fib h 2 constraint; one caveat is that systematic errors such as those due to continuum-fitting might 
have been underestimated - this is an important question which deserves closer study. In addition to better 
measurements of Jhi, future improvements in constraints on f^/i 2 from the Lya forest will rely on more 
accurate knowledge of the thermal state of the IGM. 



3. The Ionizing Continuum of Lyman Break Galaxies 

The most stringent lower limit on VL^h 2 (dashed curve in Fig. |l^i) depends sensitively on the ionizing 
continuum (E > 13.6eV) emitted by high-redshift galaxies. This is obtained from the observed composite 
spectrum of 29 Lyman break galaxies at z ~ 3.4 (Steidel et al. 2001). Because of the importance of 
the ionizing continuum in determining the likelihood for fi^/i 2 , we discuss here the results of Steidel et 
al. (2001) in some more detail. In particular, the continuum flux in the Steidel et al. (2001) sample is 
unexpectedly high (compared with both naive theoretical expectations, and observations of galaxies at low 
redshifts; see below), and, taken at face value, presents two distinct puzzles: (i) a significant fraction of the 
ionizing radiation appears to escape from the LBGs, whereas a small escape fraction is expected; and (ii) 
even assuming that no ionizing photons were lost in escaping from the galaxies, the stellar populations of 
the LBGs appear to produce a very hard ionizing continuum compared to spectral synthesis models. We 
address these two issues below, first (ii), and then (i). 



3.1. The Ionizing Flux Produced in Lyman Break Galaxies 

We first compare the stacked composite spectrum of LBGs in Steidel et al. (2001) with model galactic 
spectra. We utilize the spectral synthesis models of Leitherer et al. (1999); although we obtained similar 
conclusions using the models of Bruzual & Chariot (1996). We adopt the stellar models for a Salpeter IMF 
(with a slope a = 2.35) between 1 — 100M Q , metallicity of Z = 0.4Z Q , and compute the spectrum either 
for continuous star formation for 10 8 yr, or at the age of 10 6 yr following an instantaneous starburst. These 
correspond to the models shown in Figures 8d and 7d of Leitherer et al. (1999), respectively. In addition, to 
simulate the observed composite spectrum, we modify the emitted stellar template spectrum by including 
the effects of dust absorption, as well as the opacity of the intervening IGM. To include dust opacity, we 
use the absorption cross-sections for a mix of graphite and silicate grains that reproduces the Milky Way 
opacities (Draine & Lee 1984), and adjust the total amount of dust by fitting the slope of the observed 
spectrum at (emitted) wavelengths A > 1215 A. 

In order to assess the intrinsic level of the flux below 912A, it is important to model the opacity of the 
intervening IGM accurately. This opacity is due to Lyman a absorption systems along the lines of sights to 
the galaxies. We model here the IGM absorption by summing the optical depths of all Lyman a absorbers in 



'The Steidel et al. measurements strictly speaking constrain Jhi at z ~ 3.4 rather than 3, as we have used in §2.2 



Within 

the present uncertainties, Jjji does not appear to evolve significantly with redshift around z ~ 3 (see e.g. Giallongo et al. 1996, 
Cooke et al. 1997, Scott et al. 2000). We choose to focus on z = 3 partly because some of the best existing measurements on 
mean decrement and transmission power spectrum are at 2 = 3. 
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the range of neutral HI column densities 10 12 — 10 20 cm -2 (see Madau 1995 for a description of the method, 
and Fardal, Giroux & Shull 1998 for a more recent parameterization of the column density distribution of 
absorbers). Note that the opacity at wavelengths shorter than 912A is dominated by the poorly-known 
abundance of high-column density absorbers (Ahi 10 17 cm -2 ). However, the relevant integrated quantity 
that determines the break in the spectrum at 912A is better constrained. Using a quasar sample, Steidel et 
al. (2001) measured the value of the intrinsic ratio of the decrements at 900A and llOOA, which they find 
to be a factor of 2.5. In order to match this decrement ratio, and to simultaneously predict a Lya decrement 
of Da ~ 0.45 (which we obtained directly from the z = 3.4 LBG sample), we found that we had to adopt a 
column-density distribution intermediate between models Al and A4 of Fardal et al. (1998). 

In Fig. H| we show the resulting model spectra, superimposed on the composite spectrum of Steidel et al. 
(2001). The dotted and dashed curves correspond to the models of continuous star formation for 10 8 years, 
and to the 10 6 yr old starburst, respectively; the solid curves show the data. All fluxes are normalized to 
the same flux at the emitted wavelength of 1500A. In the bottom panel, we also show the optical depths we 
assumed for the dust and the IGM, respectively (in the continuous star formation case; a somewhat larger 
dust opacity was needed for the instantaneous starburst). In the top panel, the horizontal dashed line shows 
the mean observed flux in the wavelength range 880A < A < 910A. 

As the top panel in Fig. || shows, in the 880A < A < 910A range, the mean observed flux is F v /Fi5oo — 
0.06 ± 0.01; while the continuous star formation models predict i^/i^isoo = 0.027 ± 0.003. Taken at face 
value, the data reveals twice as high an ionizing flux as these models predict, even if we do not include any 
absorption by cold gas in the galaxy. The starburst model can come close to producing the observed ionizing 
flux at 912A, but only if an age of ^ 10 6 yr is assumed. Although the 29 galaxies included in the composite 
spectrum were selected from the bluest quartile of the sample, it appears unlikely that they were typically 
"caught" at such a young age. A short lifetime does appear to be consistent with the number density and 
clustering properties of LBG's (e.g. Wechsler et al. 2001; but see also Mo, Mao & White 1999), but it is 
disfavored by other arguments, such as the optical to near- infrared colors (Shapley et al., in preparation), 
and by tentative measurements of mass-to-light ratios (Pettini et al. 2001). The ionizing radiation in the 
models is dominated by massive OB stars, and, as a result, after ~ 4x 10 6 years, the starburst model predicts 
a rapidly diminishing value for the flux below 912A. 

The starburst model gives a (cross-section weighted) spectral slope of (3 — 1.4 below 912A (j„ cx 
This slope is somewhat harder than typical quasar spectra in the relevant wavelength range just below 912A. 



This is what we have used in §2.2 to convert measurements of j„ HI to Jm. The true value for (3 is likely to be 
smaller due to absorption by the intergalactic medium - this gives us a lower limit on Jm (Jm oc l/(/3 + 3); 
see eq. [§). 

Taking the data at face value, the implication is that the stellar populations in LBGs produce ^> twice 
as many H-ionizing photons as the current models predict. This discrepancy is more severe at shorter 
wavelengths, where the IGM is expected to become increasingly more opaque. Note that although the mean 
IGM opacity is directly measured by Steidel et al. at A « 900A, the variation of the opacity with wavelength 
across the range 880 — 910A depends on the column density distribution of Lyman a absorption systems 
with column densities iVm ^5 10 17 cm -2 , and is more uncertain. 

Changing the assumed metallicities, or adding the emission from the nebular continuum, does not 
significantly increase the predicted ionizing fluxes. The UV fluxes of OB stars can be increased by the 
presence of stellar winds (not included in the stellar models we adopted). Although winds can increase the 
Hell-ionizing flux by orders of magnitude, the corresponding increase for the H-ionizing flux in O stars 
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has been found to be small (Schaerer & de Koter 1997). The increase can be more significant for cooler 
B-stars (Najarro et al. 1996), but these stars do not dominate the ionizing photon budget in the continuous 
star-formation models, in which O-stars are continuously replenished. Furthermore, the hydrostatic stellar 
models typically reproduce the properties of Galactic HII regions (Leitherer et al. 1999), so that an increase 
by the required factor of ~ 2 would make it more difficult to reconcile the models with these observations. 

In summary, we conclude that the stellar population in the galaxies comprising the Steidel et al. (2001) 
sample appear to be significantly different from spectral synthesis model predictions, and from stellar pop- 
ulations in nearby galaxies, unless the selection procedure somehow strongly favors very young galaxies (i.e. 
age < 10 6 years). It remains to be seen whether stellar models that include winds can account for the 
apparent difference. An alternative explanation of the large ionizing flux could be that the IMF at z ~ 3 is 
biased towards massive stars. 



3.2. The Escape Fraction of Ionizing Photons From Lyman Break Galaxies 

Although the nature of LBGs is not well understood, based on their observed fluxes and clustering 
properties, they appear to have luminosities similar to galaxies, and to be associated with halos of total 
mass <; 10 11 M© (in most models, see, e.g. Wechsler et al. 2001, and references therein). The escape of 
ionizing radiation from "normal" disk galaxies has been studied extensively in observations of the Milky 
Way (Reynolds et al. 1995) and other nearby galaxies (e.g. Leitherer et al. 1995; Hoopes, Walterbos & 
Rand 1999); as well as theoretically (e.g. Dove & Shull 1994; Dove, Shull & Ferrara 2000). Both theory 
and observations indicate small escape fractions that are of order ~ 10%. In a theoretical study extending 
galactic disk models to high redshift, Wood & Loeb (2000) showed that if a few percent of the total gas 
mass would cool and settle into a rotationally supported thin disk, the escape fraction would be negligible 
from galaxies at z ~ 3 in halos with a total mass M <; 1O U M0 (the decrease in / osc being due to the higher 
densities [oc (1 + z) 3 ] at higher redshifts). In comparison, the escape fraction in the Steidel et al. sample is 

/esc £ 0.5. 

Ionizing radiation could escape more efficiently from LBGs if most of the cold gas (T ~ 10 4 K) were 
spread over an extended region, rather than incorporated into a fully-formed, dense galactic disk. Such 
a spatially extended distribution for the cold gas would be natural during the earliest phases of galaxy 
formation, when galaxy-sized dark halos are assembling for the first time. Alternatively, if star formation 
in LBGs is triggered by mergers (see this model in Wechsler et al. 2000), then the merger could spread the 
cold gas over a large solid angle. In either case, the diluted density would reduce the recombination rate 
(oc p 2 ), making it easier to keep the cold gas photoionized, and for the ionizing radiation to escape. 

To quantify the required ionizing photon rate, we consider the simplified toy model for the cold gas 
distribution described in Haiman & Rees (2001). In this model, the gas initially has a spherically symmetric 
distribution, with a radial profile dictated by hydrostatic equilibrium in an NFW (Navarro, Frenk & White 
1997) dark matter halo. Subsequently, a fraction of the gas, determined by the cooling efficiency, cools to 
10 4 K, and resides in cold clumps whose densities are compressed by a factor ~ T v i r /10 4 K (where T v „ ~ 10 6 K 
is the virial temperature of the halo). The rest of the gas remains in the collisionally ionized hot phase, at 



9 Steidel et al. define / eS c as the fraction of 900A photons that escape the galaxy, divided by the same fraction for 1500A 
photons. Thus, the actual escape fraction of 900A photons from the LBGs can be as small as ~ 15%, if dust absorbs ~ 85% 
of photons at both wavelengths. Since dust opacity is not included in the theoretical models, this is a fair comparison. 
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the temperature of T vlr . 

As an example of this model, the gas in a 10 12 M Q halo at z = 3.4 (virial temperature of T vir = 2 x 10 6 K, 
a 2.3cr object) cools out to approximately half of the virial radius (i? v ; r ~ 50kpc). The total recombination 
rate within the halo is R — Air J^ v " r 2 drfvnn(r) 2 aB is~ 2x 10 55 s _1 , where /y is the volume filling fraction 
of cold gas (approximately unity out to i? v ir/2, and dropping to negligible values at larger radii), and as is 
the case-B hydrogen recombination coefficient evaluated at « 10 4 K. With the standard Salpeter IMF, every 
stellar proton yields ~ 4000 ionizing photons, so that a star formation rate of 1 M Q yr _1 translates to a 
production rate 1.5 x 10 53 s _1 of ionizing photons. In order to allow a large escape fraction, most of this 
cold gas has to be kept photo-ionized. Increasing the above photon production rate by a factor of four, 
as suggested by the composite spectrum of Steidel et al. (2001), we find that the required "characteristic" 
star formation rate is ~ 40M Q yr _1 . This star formation rate is consistent with the luminosities of LBGs. 
Furthermore, this rate converts only ~ 2% of the available gas into stars in 10 8 yr, making our assumption 
of continuous star formation for ~ 10 s years plausible. We also note that in this model, the required star 
formation rate for a halo of mass M at redshift z scales approximately as 

making it easier for ionizing radiation to escape from smaller halos (see Haiman & Rees 2001 for details). 

The model described above is highly idealized, and should be regarded only as a plausibility argument 
that large escape fractions are possible with an extended gas distribution. These models predict that the 
photoionized cold gas could be detectable in the Lya recombination line, as a low surface brightness "fuzz" 
with <; 10 -18 erg s _1 cm~ 2 asec~ 2 . Deep Lya imaging of extended regions around the LBGs can therefore 
furnish a method to test this model (see Haiman & Rees 2001 for a similar discussion in the context of bright 
quasars). 

An alternative scenario is one in which the cold gas resides in clumps with a small covering factor. This 
could allow most of the ionizing photons to escape along lines of sights traversing only hot, collisionally 
ionized, medium. However, in a self-gravitating two-phase medium, this explanation requires a minimum 
cold clump size that exceeds the Jeans mass of the cold phase (Rees 1988). As a result, the postulated 
large clumps would be unstable to fragmentation, which would tend to increase the covering factor to unity. 
Nevertheless, the covering factor of cold gas can be decreased by feedback mechanisms. Indeed, Pettini et 
al. (2001) found evidence for galactic scale outflows ( "superwinds" of several hundred km/s) in the spectra 
of LBGs to be ubiquitous. As noted in Pettini et al. (2001), such superwinds, driven by mechanical energy 
from supernova and stellar winds, can punch holes through the galactic disk, through which the ionizing 
radiation can escape. 



4. Discussion 

To summarize, the accuracy of the present constraints on the baryon fraction from the Lya forest is 
limited by our knowledge of the ionizing background. A slight tension does exist between such constraints and 
the nucleosynthesis value, in that the forest measurements tend to prefer higher values for fl^h 2 . However, 
even assuming the largest estimate of the ionizing background, which is dominated by contributions from 
Lyman-break galaxies as reported by Steidel et al. (2001) (see model 3 in § 2.2.1 ), the nucleosynthesis value 
is ruled out only at the 3a level. 
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We emphasize also that even if the ionizing background is known to high accuracy, substantial uncer- 
tainties remain, largely due to the uncertain thermal state of the intergalactic medium (parametrized by Tq 



and 7 - the temperature at mean density and the equation of state index; see §2.1). Present measurements 
give a probability distribution of Q, b h 2 which is well-approximated by a Gaussian: 

P*(n b h 2 ) = —=J- exp [-(n b h 2 - 0.046) 2 /0.0061 2 ] (8) 

V27r0.0061 



if Jhi is fixed to be unity. We show in eq. (^) how the above probability can be used to obtain the probability 
distribution for fl b h 2 , P(flbh 2 ), given any distribution for Jhi (in the extreme case, a delta function). This 
will be useful as measurements of Jhi improve. We also show in §|2.2.1 how the lower bound likelihood can 



be obtained from a given P(flbh 2 ) - this is a more appropriate quantity to consider if only strictly lower 
limits to the ionizing background are known. 

Our results show that future improvements in the Lya forest constraint on the baryon fraction will have 
to rely on improvements in measurements of the ionizing background and of the thermal state of the gas. 
The former can be achieved through a variety of means. Measurements of the ionizing flux from LBG's 
extending to the redder population and fainter magnitudes will be useful. Improved measurement of the 
proximity effect would, in principle, be especially powerful because it measures the sum total background of 
all sources. However, one main source of uncertainty is the large velocity shifts between different lines which 
affects estimates of the effect as a function of distance from quasars (see e.g. Scott et al. 2000). In addition, 
systematic effects due to clustering and peculiar velocities should be taken into account in the translation 
of absorption line statistics into the background flux - interestingly, Loeb & Eisenstein (1995) argued that 
these effects generally lead to an over-estimate of Jhi from standard proximity effect measurements. It is also 
useful to keep in mind the possibility of measuring Jhi from direct imaging of the fluorescent Lya photons 
from optically thick systems (Gould & Weinberg 1996). The measurements of Bunker, Marleau & Graham 
(1998) have already put an upper limit on the ionizing background, which at the moment is still quite large 
(3-sigma upper limit is jhi < 2 x 10~ 21 ergs -1 cm -2 Hz -1 sr _1 ). Deeper exposure to seek a definite detection 
would be very useful. As far as the thermal state of the intergalactic medium is concerned, larger sample 
size is crucial, and several techniques have proven useful for this purpose (e.g. Ricotti, Gnedin & Shull 2000, 
Bryan et al. 1999, Schaye et al. 1999, McDonald et al. 2001, Choudhury et al. 2000, Zaldarriaga et al. 2000; 
but see also Meiksin, Bryan & Machacek 2001). We should also mention that a constraint on the baryon 
density can be obtained from the He II forest as well (e.g. Wadsley, Hogan & Anderson 2000). 

Ultimately, as the error-bars on flbh 2 from the forest come down, we will have a very precise test of the 
gravitational instability model, for the forest in particular as well as for large scale structure in general. A 
statistically significant disagreement among constraints from the forest, nucleosynthesis and the microwave 
background would require some basic revision of our model. 

The findings of Steidel et al. (2001) regarding the LBG Lyman-continuum deserve further study. Their 
measurements imply an intrinsically hard stellar spectrum, which seems to require a (perhaps unreasonably) 
young population of stars (1 million years or so) and/or an initial mass function that is biased towards 
massive stars. Moreover, the large escape fraction for ionizing photons imply that the cold, optically thick 
gas in LBG's is likely distributed quite differently from what has been seen in local galaxies. Possibilities 
include a model in which the cold gas is spatially extended over the halo, and has not settled into a dense 
galactic disk; or a disk with large holes punched out by strong stellar winds, as suggested by Pettini et al. 
(2001). 

Finally, it is interesting to note that there are some indications of a similar evolutionary trend in AGN 
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- a decrease of the escape fraction at lower redshifts. Searches for intrinsic absorption at the Lyman-edge 
at intermediate redshifts (0.4 <, z <J 2) generally turn up only small Lyman edge discontinuities in a small 
fraction of quasars (see e.g. Koratkar, Kinney & Bohlin 1992, Zheng et al. 1997). On the other hand, there 
is mounting evidence of strong internal Lyman edge absorption in nearby Seyfert galaxies (Alexander et al. 
1999, 2000; Kraemer et al. 1999). 

We thank C. Leitherer and M. Pettini for useful discussions, and C. Steidel for electronic data for the 
observational spectrum in Fig. |^. LH was supported by an Outstanding Junior Investigator Award from 
DOE. ZH was supported by NASA through the Hubble Fellowship grant HF-01119.01-99A, awarded by 
the Space Telescope Science Institute, which is operated by the Association of Universities for Research in 
Astronomy, Inc., for NASA under contract NAS 5-26555. 
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Fig. 1. — (a) Lower panel: likelihood of A(To/10 4 K) ' at z = 3. (b) Upper panel: lower limit to the 
cumulative probability that VL^h 2 is larger than some value Bq (eq. Dotted, solid and dashed lines 

correspond to choices 1, 2 and 3 for P( Jhi) respectively (see text). The horizontal bar labeled by BBN gives 
the nucleosynthesis 68% region. 
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Fig. 2. — (a) Lower panel: the likelihood (differential rather than cumulative distribution) for flbh 2 , where 
dotted, solid and dashed lines correspond to choices 1, 2 and 3 for P(Jhi) respectively. These curves derive 
from differentiation of those in Fig. [j]b; see text for details, (b) Upper panel: P*(f2;,/i 2 ), the likelihood for 
f^ft 2 if Jjji is fixed to be 1; this is useful for generating the likelihood of Sl^ti 2 for arbitrary distributions of 
Jhi (sec cq. ||). P^Vt^h 2 ) is well approximated by a Gaussian of mean = 0.046 and dispersion = 0.0061. 
For Jhi fixed to be some value different from unity, the probability distribution for f^/i 2 can be obtained by 
simply multiplying both the mean and dispersion by \fJm_. 
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Fig. 3. — Population synthesis spectra, including opacities of dust and neutral hydrogen in the IGM, 
compared with the observed spectrum of Steidel et al. (2001). The dotted and dashed curves correspond to 
models of continuous star formation for 10 s years, and to an instantaneous starburst at 10 6 yr, respectively. 
The solid curves show the Steidel et al. (2001) composite spectrum. All fluxes are normalized to the flux at 
the emitted wavelength of 1500A. For reference, in the bottom panel, we show the optical depths assumed 
for the dust and the IGM. In the top panel, the horizontal dashed line shows the mean observed flux in the 
interval 880l < A < 910l. 



